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Abstract 

A model equation for the motion of a vortex filament immersed in three dimensional, 
incompressible and inviscid fluid is investigated as a humble attempt to model the motion 
of a tornado. We solve an initial-boundary value problem in the half space where we impose 
^ ■ a boundary condition in which the vortex filament is allowed to move on the boundary. 

1 Introduction 

Many researchers have studied tornadoes from several perspectives. A systematic research and 
observation of tornadoes is difficult mainly because of two reasons: a precise prediction of 
tornado formation is not yet possible, and the life-span of a tornado is very short, giving only 
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short openings for any kind of measurements. Many aspects of tornadoes are still unknown. 

In 1971, Fujita [3] gave a systematic categorization of tornadoes. He proposed the so-called 
Fujita scale in which tornadoes are classified according to the damage that it dealt to buildings 
and other surroundings. The scale provides a correlation between ranges of wind speed and 
the damage that it causes. The enhanced version, called the Enhanced Fujita Scale, is used to 
classify tornadoes to date. McDonald |10j gives a review of Fujita's contributions to tornado 
research. 

Since then, due to the advancement of technology, more accurate and thorough observations 
and simulations have become possible, and theories for the formation and motion of tornadoes 
have developed. Klemp [5] and the references within give an extensive review on the known 
dynamics of tornadoes. 

Motivated by this, we investigate the motion of a vortex filament. The vortex filament 
equation, also called the Localized Induction Equation (LIE) models the movement of a vortex 
filament, which is a space curve where the vorticity of the fluid is concentrated, and is described 
by 

(1.1) = X s X X SS) 

where x(s,t) = (x 1 (s,t), x 2 (s,t), x 3 (s,t)) is the position of the vortex filament parameterized 
by the arc length s at time t, x denotes the exterior product, and the subscripts denote differ- 
entiation with respect to that variable. We also use d s and dt for partial differentiation with the 
corresponding variables. 

The LIE was first derived by Da Rios [2] and re-examined by Arms and Hama pQ. Since 
then, many authors have worked with the equation. Nishiyama and Tani 18| gave the unique 
solvability of the initial value problem for the LIE in Sobolev spaces. A different approach was 
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taken by Hasimoto [3J. He used the so-called Hasimoto transformation to transform (jl.ip into 
a nonlinear Schrodinger equation: 

ld^ d*i> 1 2 

where ^ is given by 

ip = k exp I i / r ds 



k is the curvature, and r is the torsion of the filament. Even though this expression is undefined 
at points of the filament where the curvature vanishes, Koiso [6] proved that the Hasimoto 
transformation is well-defined in the class of C°° functions. He used a geometrical approach 
to define the Hasimoto transformation and showed the unique solvability of the initial value 
problem in the class of C°° functions. 

Regarding initial-boundary value problems, the only known result that the authors know is 
by Nishiyama and Tani [8 J . The boundary condition imposed there necessarily fixes the end point 
of the vortex filament and does not allow it to move on the boundary. From the physical point of 
view, the vortex filament must be closed, extend to the spatial infinity, or end on boundaries of 
the fluid region. In the last case, we have to impose an appropriate boundary condition to show 
the well-posedness of the problem. Since it is hard to find what kind of boundary condition is 
physically reasonable if we begin our analysis from the Schrodinger equation, we chose to work 
with the original vortex filament equation (jl.ip . 

In light of modeling the motion of a tornado, we consider (jl.ip in a framework in which the 
end of the vortex filament is allowed to move on the boundary. We do this by setting a different 
boundary condition than that of [S|. 

The contents of this paper are as follows. In section 2, we formulate our problem and give 
basic notations. In section 3, we derive compatibility conditions for our initial-boundary value 
problem. Sections 4 and 5 are concerned with constructing the solution. 



2 Setting of the Problem 

We consider the initial-boundary value problem for the motion of a vortex filament in the half- 
space in which the filament is allowed to move on the boundary: 

{x t = x s x x ss , s > 0, t > 0, 
x{s,0) = x {s), s>0, 
x s (0,t)=e 3 , t>0, 

where e% = (0, 0, 1). We assume that 

(2.2) |£c 0s (s)| = 1 for s > 0, a;jj(0) = 0, 

for the initial datum. The first condition states that the initial vortex filament is parametrized 
by the arc length and the second condition just states that the curve is parameterized starting 
from the boundary. Here we observe that by taking the inner product of e 3 with the equation, 
taking the trace at s = 0, and noting the boundary condition we have 

jjj (e 3 • x) \ s=0 = e 3 • (x s x x ss )\ s=Q 

= X s • {x s X 33 ss )| s _q 

= 0, 
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where "■" denotes the inner product and | s= o denotes the trace at s = 0. This means that if the 
end of the vortex filament is on the boundary initially, then it will stay on the boundary, but is 
not necessarily fixed. This is our reason for the notion "allowed to move on the boundary" . 
By introducing new variables v(s,t) := x s (s,t) and vq(s) := xq s (s), (|2.ip and (|2.2p become 

{v t = v x v ss , s > 0, t > 0, 

t>(s,0) =«„(*), s>0, 
v(0,t) = e 3 , t>0, 



(2.4) Ms)| = l, s>0. 

Once we solve (|2.3p . the solution x of (|2.ip and (|2.2p can be constructed by 

x(s,t) = x (s) + / v(s,t) X V s (s,T)dT. 
Jo 

So from now on, we concentrate on the initial-boundary value problem (|2.3|) under the condition 
(|2.4p . Note that if the initial datum satisfies (j2.4j) , then any smooth solution v of (|2.3p satisfies 

(2.5) KM)I = 1, a > 0, t > 0. 

This can be confirmed by taking the inner product of the equation with v. 
We define basic notations that we will use throughout this paper. 

For a domain f2, a non-negative integer m, and 1 < p < oo, W m ' p (Q) is the Sobolev space 
containing all real- valued functions that have derivatives in the sense of distribution up to order 
m belonging to L p (f2). We set H m (Q) = W m ' 2 (Q) as the Sobolev space equipped with the usual 
inner product, and H°(Q) = L 2 (Q). We will particularly use the cases Q = R and Q = R+, 
where R+ = {s G R; s > 0}. The norm in H m (Q) is denoted by || • || m and we simply write || • || 
for || • 1 1 o - We do not indicate the domain in the symbol for the norms since we use it in a way 
where there is no risk of confusion. 

For a Banach space X, C m ([0,T]; X) denotes the spaces of functions that are m times 
continuously differentiable in t with respect to the topology of X. 

For any function space described above, we say that a vector valued function belongs to the 
function space if each of its components does. 



3 Compatibility Conditions 

We derive necessary conditions for a smooth solution to exist for (|2.3p with (12. 4p . 

Suppose that v(s,t) is a smooth solution of (|2.3p with (|2,4p defined in R + x [0, T] for some 
positive T. We have already seen that for all (s,t) G R+ x [0, T] 

(3.1) \v(s,t)\ 2 = l. 

By differentiating the boundary condition with respect to t we see that 

(.B)„ ofu| s=0 = for nGN, t > 0. 

We next show 
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Lemma 3.1 For a smooth solution v(s,t) under consideration, it holds that 



(C) n 

(D) n 



v x d 2n v 



0, 



d J s v ■ d s v 



s=0 



s=0 

for j + l = 2n + l. 



Proof. We prove them by induction. From (-B)i and by taking the trace of the equation we see 
that 

= v t \ s=0 = vxv ss \ s=0 , 
thus, (C)i holds. By taking the exterior product of v s and (C)i we have 

{(v s ■ v ss ) v - (v s ■ v) v ss } \ s=0 = 0. 

On the other hand, by differentiating (13. ip with respect to s we have v ■ v s = 0. Combining 
these two and the fact that v is a non-zero vector, we arrive at 

v s ■ v ss \ s=0 = 0. 

Finally, by differentiating (|3.1|) with respect to s three times and setting s = 0, we have 

= 2 (v ■ v sss + 3v s ■ v ss ) \ S=Q = 2v ■ v sss \ s=0 , 

so, (-D)i holds. 

Suppose that the statements hold up to n — 1 for some n > 2. By differentiating (C) n _i 
with respect to t we have 

»x(^("- 1 )t; t )| a=0 = J 
where we have used (-B)i- We see that 

2(n-l) 



where 



(3.2) 



d^v t = ^""D (v x v ss ) = ( 2(n fc " 1} ) W« x d 2 s ^- k+2 v), 

2(n - 1) 



is the binomial coefficient. So we have 



2(n-l) 

£ 

A:=0 



2(n-i; 
k 



0. 



s=0 



We examine each term in the summation. When 2 < k < 2(n — 1) is even, we see from the 
assumptions of induction (C%/2 an d (C)(2(n-i)-fc+2)/2 that both d k v and <9 2 ^ n k+2 v are 
parallel to u, so that 

d k s v x a s 2(n - 1} - fc+2 t; = 0. 

s=Q 

When 1 < k < 2(n — 1) is odd, we rewrite the exterior product in (|3.2p as 

« x (d k v x 9 s 2 ("- 1 )- fc + 2 - u ) = (« • df n -^- k+2 v)d k v - (v ■ d k v)df n -^- k+2 v. 

Since 2(n — 1) — + 2 is also odd, by (-D)(fc-i)/2 and (-D)(2(n-i)-fc+i)/2 we have 



v ■ d 2 S n -^- k+2 v 



s=0 



0. 



s=0 
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Thus, only the term with k = remains and we get 

v x (v x d 2n v)\ s=Q = 0. 

Here, we note that 



where we used (|3.ip . Taking the exterior product of this with v we see that (C) n holds. Taking 
the exterior product of d 2n+1 ~ 2k v w ^ (C)k and using {D) n _k for 1 < k < n yields 



(dfv ■ d 2n+l ~ 2k v)v 
Since v is non zero, we have for 1 < k < n 



= 0. 

s=0 



(3.3) dfv ■ d 2 s n+1 - 2k v 



= 0. 

s=0 



Finally, by differentiating (|3.ip with respect to s (2n + 1) times, we have 

2n+l 



E ( / ) ( d > ■ % n+1 ~ jv ) 



= 0. 

s=0 



3=0 

Since every term except j = 0, 2n + 1 is of the form (|3.3p . we see that 

vd 2n+1 v\ =0, 

* \s=0 ' 

which, together with (|3.3p . finishes the proof of (.D) n . □ 

Worth noting are the following two properties which will be used in later parts of this paper. 
For integers n, 

e 3 x d 2n v\ =0, e 3 ■ 5 s 2n+1 ^| _ = 0. 

These are special cases of (C) n and (-D) n with the boundary condition substituted in. 

By taking the limit t — > in (C) n , we derive a necessary condition for the initial datum. 

Definition 3.2 For n £ N U {0}, we say that the initial datum vq satisfies the compatibility 
condition (A) n if the following condition is satisfied for < k < n 

j v o\ s=0 = e 3 , k = 0, 

\ {v xd 2k v )\ s=Q = 0, fc£N. 

From the proof of Lemma [3. H we see that if vq satisfies (|2.4p and the compatibility condition 
(A) n , then vq also satisfies (-D)fc for < k < n with v replaced by vq as long as the trace exists. 

4 Extension of the Initial Datum 

For the initial datum v$ defined on the half-line, we extend it to the whole line by 



(4.1) v (s) = 

where v = {v l ,v 2 , — u 3 ) for v = (v , v 2 , v 3 ) € R 3 . 



v (s), s>0, 
-vo(-s), s < 0, 
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-(-l) fc+1 (d s fc+ V)(- S ), s<0, 



Proposition 4.1 For any integer m > 2, ifvQ S G H m (R + ) satisfies (|2.4p and the compatibility 
condition (A)[m], i/ien ^os G H m (R.). Here, [y] indicates the largest integer not exceeding y. 

Proof. Fix an arbitrary integer m > 2. We will prove by induction on k that 8^ +1 vq G £ 2 (R) 
for any < k < m. Specifically we show that the derivatives of Vq in the distribution sense on 
the whole line R up to order m + 1 have the form 

f (d* +1 v ) (a), s > 0, 

(4.2) (d k s +1 v )(s)= l [S ° >{h 

for < k < m. 

Since vq G L°°(R + ) and i?o s G i/ 2 (R + ), Sobolev's embedding theorem states t^os G L°°(R + ) 
and thus vq G Ty 1,00 (R+), so that the trace vq(0) exists. By definition (|4.ip we have 

«o(-0) = (-«5(0), - Vo 2 (0), « 3 (0)) , 

but from (-A)o, ^o(O) = Vq(0) = 0. These imply that vq(+0) = vq(—0), so that we obtain 

/ (d s v ) (a), s > 0, 

\ -(-l)(i\^)(-s), s<0, 



d s v (s) 



and the case k = is proved. 

Suppose that (I4.2p with k + 1 replaced by k holds for some k £ {1,2, ... , m}. We check that 
the derivative O^vq does not have a jump discontinuity at s = 0. When is even, from the 
definition of O^vq, 

(#t>„)(-0) = (-«(0), -0*ug(O), ^ 3 (0)), 

but from we have 



■i 



= v x <9> = e 3 x 3> (0), 

s=0 

which means that dgVo(0) is parallel to e 3 and that the first and second components are zero. 
When k is odd, 

(flS)(-o) = (fl*«g(o), a^ 2 (o), -a^ 3 (o)), 

but (-A)rfci implies (-D)rfci and particularly 



= v ■ d k s v Q 



s=0 



e 3 -5%(0) = (^ 3 )(0), 



so the third component is zero. In both cases, we have (dgVo)(+Q) = (dgVtj) (— 0), so that we 
can verify ()4.2p . This finishes the proof of the proposition. □ 

5 Existence and Uniqueness of Solution 

Using So j we consider the following initial value problem: 

(5.1) u t = u x u ss , s G R, t > 0, 

(5.2) u(s,0) = v (s), sGR. 

By Proposition ^. 11 the existence and uniqueness theorem (cf. Nishiyama [9]) of a strong solution 
u is applicable. Specifically we use the following theorem. 
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Theorem 5.1 (Nishiyama [9]) For a non-negative integer m, if vq s E H 2+m (R) and \vq\ = 1, 
then the initial value problem (|5.ip and (|5.2p has a unique solution u such that 



u-v £ C([0,oo);^ 3+m (R)) nC 1 ([0,oo);i7 1+m (R)) 

and \u\ = 1. 

From Proposition 14.11 the assumptions of the theorem are satisfied if vq s E H 2+m (R + ) 
satisfies the compatibility condition (A) ] 2+ m ^ and (|2.4p . 
Now we define the operator T by 

WW = -»(-*), 

for R 3 -valued functions io defined on s E R. By direct calculation, we can verify that Ti;o = ?>o 
and that T (it x tt ss ) = (Tit) x (Tu) ss . Taking these into account and applying the operator T 
to ([53]) and (IQ]) . we have 

f (Ttt)t = (T«) x (Tt*) M> s E R, i > 0, 

| (Ttt)(a, 0) = (Tt>o)(s) = »o(«), a G R, 

which means that Tu is also a solution of (|5.ip and (|5.2p . Thus we have Tu = it by the 
uniqueness of the solution. Therefore, for any t E [0, T] 

u(0,t) = (Tu) (0,t) = -u{0,t), 

which is equivalent to u 1 (0, t) = u 2 (0, t) = 0. Therefore, it holds that u 3 (0, i) = — 1 or 1 because 
we have |u| = 1. But in view of «o(0) = fo(0) = e3, we obtain u(0, t) = e% by the continuity in 
t. 

This shows that the restriction of u to R+ is a solution of our initial-boundary value problem. 
Using this function v := w|r + , we can construct the solution x to the original equation as we 
stated in section 2. Thus we have 

Theorem 5.2 For a non-negative integer in, if xq ss E H 2+m (R + ) and xq s satisfies the com- 
patibility condition (A)^ 2+ m ^ and ()2.2fl . then there exists a unique solution x of (12. ip such that 

x-x e C([0,oo);F 4+m (R + )) nC 1 ([0,oo); J ff 2+m (R+)), 

and \x s \ = 1. 

Proof. The uniqueness is left to be proved. Suppose that X\ and xi are solutions as in the 
theorem. Then, by extending ajj (i = 1,2) by 



£j(s,i) 



Xi(s,t) s > 0, i > 0, 
Xi(-s,t) s < 0, i > 0, 



we see that ajj are solutions of the vortex filament equation in the whole space. Thus x\ = x<i 
follows from the uniqueness of the solution in the whole space. □ 
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